Simultaneous species detection and discovery with environmental DNA metabarcoding: A freshwater mollusk case study.

Environmental DNA (eDNA) sampling is a powerful tool for rapidly characterizing biodiversity patterns for specious, cryptic taxa with incomplete taxonomies. One such group that are also of high conservation concern are North American freshwater gastropods. In particular, springsnails of the genus Pyrgulopsis (Family: Hydrobiidae) are prevalent throughout the western United States where >140 species have been described. Many of the described species are narrow endemics known from a single spring or locality, and it is believed that there are likely many additional species which have yet to be described. The distribution of these species across the landscape is of interest because habitat loss and degradation, climate change, groundwater mining, and pollution have resulted in springsnail imperilment rates as high as 92%. Determining distributions with conventional sampling methods is limited by the fact that these snails are often <5 mm in length with few distinguishing morphological characters, making them both difficult to detect and to identify. We developed an eDNA metabarcoding protocol that is both inexpensive and capable of rapid, accurate detection of all known Pyrgulopsis species. When compared with conventional collection techniques, our pipeline consistently resulted in detection at sites previously known to contain Pyrgulopsis springsnails and at a cost per site that is likely to be substantially less than the conventional sampling and individual barcoding that has been done historically. Additionally, because our method uses eDNA extracted from filtered water, it is non-destructive and suitable for the detection of endangered species where "no take" restrictions may be in effect. This effort represents both a tool which is immediately applicable to taxa of high conservation concern across western North America and a case study in the broader application of eDNA sampling for landscape assessments of cryptic taxa of conservation concern.

The unknown unknown: A framework for assessing environmental DNA assay specificity against unsampled taxa.

Taxon-specific quantitative PCR (qPCR) assays are commonly used for environmental DNA sampling-based inference of animal presence. These assays require thorough validation to ensure that amplification truly indicates detection of the target taxon, but a thorough validation is difficult when there are potentially many non-target taxa, some of which may have incomplete taxonomies. Here, we use a previously published, quantitative model of cross-amplification risk to describe a framework for assessing qPCR assay specificity when there is missing information and it is not possible to assess assay specificity for each individual non-target confamilial. In this framework, we predict assay specificity against unsampled taxa (non-target taxa without sequence data available) using the sequence information that is available for other confamilials. We demonstrate this framework using four case study assays for: (1) An endemic, freshwater arthropod (meltwater stonefly; Lednia tumana), (2) a globally distributed, marine ascidian (Didemnum perlucidum), (3) a continentally distributed freshwater crustacean (virile crayfish; Faxonius virilis, deanae and nais species complex) and (4) a globally distributed freshwater teleost (common carp; Cyprinus carpio and its close relative C. rubrofuscus). We tested the robustness of our approach to missing information by simulating application of our framework for all possible subsamples of 20-all non-target taxa. Our results suggest that the modelling framework results in estimates which are largely concordant with observed levels of cross-amplification risk using all available sequence data, even when there are high levels of data missingness. We explore potential limitations and extensions of this approach for assessing assay specificity and provide users with an R Markdown template for generating reproducible reports to support their own assay validation efforts.

An inside "beak": Molecular analysis of swab samples reveals the seabird diet of invasive Barn Owls in Hawai'i.

Predation is an important species interaction to monitor when assessing an invasive species' impact on a particular ecosystem, but it can be difficult to observe and thus, fully understand. On Kaua'i island, invasive Barn Owls (Tyto alba) predate native seabirds, but difficult terrain in this region and the cryptic nature of owl predation make traditional monitoring of predation quite challenging. Using Barn Owls collected as part of removal efforts on Kaua'i and Lehua islands, we conducted DNA metabarcoding of owl digestive tracts to detect and determine seabird species they predate. We used a seabird-targeted 12s marker to sequence 112 swabs from 55 owls and detected six seabird species, including two ESA-listed seabirds - Hawaiian Petrel (Pterodroma sandwichensis) and Newell's Shearwater (Puffinus newelli), in 12 swabs from 11 owls (20% of sampled owls). Corresponding morphological assessment of owl stomach contents detected seabird species as prey items in only 2% (1/55) of sampled owls, highlighting the utility of molecular approaches for detecting diet items, especially degraded or visually absent items. Additionally, this approach has proven very useful in revealing cryptic trophic interactions in inaccessible seabird populations. For the most comprehensive analysis of diet, the use of both esophageal and cloacal swabs for metabarcoding is recommended. Supplementing metabarcoding with other methods that can provide complementary prey information, such as stable isotope analysis, would help to characterize trophic interactions more fully. The method described here has proven to be a reliable tool for investigating diet in invasive owls and may be used to investigate cryptic predation in living birds as a minimally invasive technique, as well.

Fire-driven animal evolution in the Pyrocene.

Fire regimes are a major agent of evolution in terrestrial animals. Changing fire regimes and the capacity for rapid evolution in wild animal populations suggests the potential for rapid, fire-driven adaptive animal evolution in the Pyrocene. Fire drives multiple modes of evolutionary change, including stabilizing, directional, disruptive, and fluctuating selection, and can strongly influence gene flow and genetic drift. Ongoing and future research in fire-driven animal evolution will benefit from further development of generalizable hypotheses, studies conducted in highly responsive taxa, and linking fire-adapted phenotypes to their underlying genetic basis. A better understanding of evolutionary responses to fire has the potential to positively influence conservation strategies that embrace evolutionary resilience to fire in the Pyrocene.

A multidisciplinary approach to estimating wolf population size for long-term conservation.

The wolf (Canis lupus) is among the most controversial of wildlife species. Abundance estimates are required to inform public debate and policy decisions, but obtaining them at biologically relevant scales is challenging. We developed a system for comprehensive population estimation across the Italian alpine region (100,000 km2 ), involving 1513 trained operators representing 160 institutions. This extensive network allowed for coordinated genetic sample collection and landscape-level spatial capture-recapture analyses that transcended administrative boundaries to produce the first estimates of key parameters for wolf population status assessment. Wolf abundance was estimated at 952 individuals (95% credible interval 816-1120) and 135 reproductive units (i.e., packs) (95% credible interval 112-165). We also estimated that mature individuals accounted for 33-45% of the entire population. The monitoring effort was spatially estimated thereby overcoming an important limitation of citizen science data. This is an important approach for promoting wolf-human coexistence based on wolf abundance monitoring and an endorsement of large-scale harmonized conservation practices.

Una estrategia multidisciplinaria para la estimación del tamaño poblacional de los lobos para la conservación a largo plazo Resumen El lobo (Canis lupus) está entre las especies de fauna más controversiales. Se requieren estimaciones de abundancia para informar al debate público y las decisiones políticas, pero es un reto obtenerlos en escalas con relevancia biológica. Desarrollamos un sistema para la estimación completa de la población en la región alpina de Italia (100,000 km2 ), con la participación de 1,513 operadores entrenados que representan a 160 instituciones. Esta red extensa permitió una colecta coordinada de muestras genéticas y análisis de captura-recaptura espacial que trascendieron las fronteras administrativas para así producir las primeras estimaciones de los parámetros clave para la evaluación del estado de la población de los lobos. Se estimó la abundancia en 952 individuos (95% intervalo de confianza 816-1120) y 135 unidades reproductivas (es decir, manadas) (95% intervalo de confianza 112-165). También estimamos que los individuos maduros representaban el 33-45% de toda la población. El esfuerzo de monitoreo se estimó espacialmente, por lo que sobrepasó una limitación importante de la ciencia ciudadana. Esta estrategia es importante para promover la coexistencia entre lobos y humanos con base en el monitoreo de la abundancia y el apoyo a las prácticas armonizadas de conservación a gran escala.

The ties that bind the sagebrush biome: integrating genetic connectivity into range-wide conservation of greater sage-grouse.

Conserving genetic connectivity is fundamental to species persistence, yet rarely is made actionable into spatial planning for imperilled species. Climate change and habitat degradation have added urgency to embrace connectivity into networks of protected areas. Our two-step process integrates a network model with a functional connectivity model, to identify population centres important to maintaining genetic connectivity then to delineate those pathways most likely to facilitate connectivity thereamong for the greater sage-grouse (Centrocercus urophasianus), a species of conservation concern ranging across eleven western US states and into two Canadian provinces. This replicable process yielded spatial action maps, able to be prioritized by importance to maintaining range-wide genetic connectivity. We used these maps to investigate the efficacy of 3.2 million ha designated as priority areas for conservation (PACs) to encompass functional connectivity. We discovered that PACs encompassed 41.1% of cumulative functional connectivity-twice the amount of connectivity as random-and disproportionately encompassed the highest-connectivity landscapes. Comparing spatial action maps to impedances to connectivity such as cultivation and woodland expansion allows both planning for future management and tracking outcomes from past efforts.

Variation in adult and pup wolf diets at natal den sites is influenced by forest composition and configuration.

Although wolves are wide-ranging generalist carnivores throughout their life cycle, during the pup-rearing season wolf activity is focused on natal den sites where pup survival depends upon pack members provisioning food. Because prey availability is influenced by habitat quality within the home range, we investigated the relative importance of prey species for adults and pups and further examined the relationship between habitat characteristics, wolf diet, and litter size on Prince of Wales Island (POW) in Southeast Alaska. During 2012-2020, we detected 13 active den sites within the home ranges of nine wolf packs. We estimated minimum pup counts using motion-detecting cameras and individual genotypes from noninvasive samples (hair: n = 322; scat: n = 227) and quantified wolf diet composition using fecal DNA metabarcoding (n = 538). We assessed habitat composition, configuration, and connectivity within denning and annual home ranges estimated using wolf GPS-collar data. Contrary to expectations, wolves had a more constricted diet during denning season (April 15-July 31), and within this season pups had a narrower dietary niche (species richness [S] = 4) focused more on deer (relative frequency of occurrence [O/I] = 0.924) than adults (S = 15; deer O/I = 0.591). Litter size had a positive relationship with the relative frequency of deer in a wolf pack's diet. Wolf consumption of deer was positively associated with the proportion of young-growth forest (≤25 years old) within denning and annual home ranges. High levels of vegetation patch interspersion, and the density of closed logging roads were also important predictors, suggesting these habitat qualities were influential for increasing the availability of deer to wolves. Our results contrast with previous research indicating wolf pup diets included more alternate prey (i.e., beaver) than adults and emphasize the importance of deer to wolf viability on POW, especially during denning season.

Broad-scale eDNA sampling for describing aquatic species distributions in running waters: Pacific lamprey Entosphenus tridentatus in the upper Snake River, USA.

One of the most fundamental yet challenging tasks for aquatic ecologists is to precisely delineate the range of species, particularly those that are broadly distributed, require specialized sampling methods, and may be simultaneously declining and increasing in different portions of their range. An exemplar is the Pacific lamprey Entosphenus tridentatus, a jawless anadromous fish of conservation concern that is actively managed in many coastal basins in western North America. To efficiently determine its distribution across the accessible 56,168 km of the upper Snake River basin in the north-western United States, we first delimited potential habitat by using predictions from a species distribution model based on conventionally collected historical data and from the distribution of a potential surrogate, Chinook salmon Oncorhynchus tshawytscha, which yielded a potential habitat network of 10,615 km. Within this area, we conducted a two-stage environmental DNA survey involving 394 new samples and 187 archived samples collected by professional biologists and citizen scientists using a single, standardized method from 2015 to 2021. We estimated that Pacific lamprey occupied 1875 km of lotic habitat in this basin, of which 1444 km may have been influenced by recent translocation efforts. Pacific lamprey DNA was consistently present throughout most river main stems, although detections became weaker or less frequent in the largest and warmest downstream channels and near their headwater extent. Pacific lamprey were detected in nearly all stocked tributaries, but there was no evidence of indigenous populations in such habitats. There was evidence of post-stocking movement because detections were 1.8-36.0 km upstream from release sites. By crafting a model-driven spatial sampling template and executing an eDNA-based sampling campaign led by professionals and volunteers, supplemented by previously collected samples, we established a benchmark for understanding the current range of Pacific lamprey across a large portion of its range in the interior Columbia River basin. This approach could be tailored to refine range estimates for other wide-ranging aquatic species of conservation concern.

New strategies for characterizing genetic structure in wide-ranging, continuously distributed species: A Greater Sage-grouse case study.

Characterizing genetic structure across a species' range is relevant for management and conservation as it can be used to define population boundaries and quantify connectivity. Wide-ranging species residing in continuously distributed habitat pose substantial challenges for the characterization of genetic structure as many analytical methods used are less effective when isolation by distance is an underlying biological pattern. Here, we illustrate strategies for overcoming these challenges using a species of significant conservation concern, the Greater Sage-grouse (Centrocercus urophasianus), providing a new method to identify centers of genetic differentiation and combining multiple methods to help inform management and conservation strategies for this and other such species. Our objectives were to (1) describe large-scale patterns of population genetic structure and gene flow and (2) to characterize genetic subpopulation centers across the range of Greater Sage-grouse. Samples from 2,134 individuals were genotyped at 15 microsatellite loci. Using standard STRUCTURE and spatial principal components analyses, we found evidence for four or six areas of large-scale genetic differentiation and, following our novel method, 12 subpopulation centers of differentiation. Gene flow was greater, and differentiation reduced in areas of contiguous habitat (eastern Montana, most of Wyoming, much of Oregon, Nevada, and parts of Idaho). As expected, areas of fragmented habitat such as in Utah (with 6 subpopulation centers) exhibited the greatest genetic differentiation and lowest effective migration. The subpopulation centers defined here could be monitored to maintain genetic diversity and connectivity with other subpopulation centers. Many areas outside subpopulation centers are contact zones where different genetic groups converge and could be priorities for maintaining overall connectivity. Our novel method and process of leveraging multiple different analyses to find common genetic patterns provides a path forward to characterizing genetic structure in wide-ranging, continuously distributed species.

eDNAssay: A machine learning tool that accurately predicts qPCR cross-amplification.

Environmental DNA (eDNA) sampling is a highly sensitive and cost-effective technique for wildlife monitoring, notably through the use of qPCR assays. However, it can be difficult to ensure assay specificity when many closely related species co-occur. In theory, specificity may be assessed in silico by determining whether assay oligonucleotides have enough base-pair mismatches with nontarget sequences to preclude amplification. However, the mismatch qualities required are poorly understood, making in silico assessments difficult and often necessitating extensive in vitro testing-typically the greatest bottleneck in assay development. Increasing the accuracy of in silico assessments would therefore streamline the assay development process. In this study, we paired 10 qPCR assays with 82 synthetic gene fragments for 530 specificity tests using SYBR Green intercalating dye (n = 262) and TaqMan hydrolysis probes (n = 268). Test results were used to train random forest classifiers to predict amplification. The primer-only model (SYBR Green results) and full-assay model (TaqMan probe-based results) were 99.6% and 100% accurate, respectively, in cross-validation. We further assessed model performance using six independent assays not used in model training. In these tests the primer-only model was 92.4% accurate (n = 119) and the full-assay model was 96.5% accurate (n = 144). The high performance achieved by these models makes it possible for eDNA practitioners to more quickly and confidently develop assays specific to the intended target. Practitioners can access the full-assay model online via eDNAssay (https://NationalGenomicsCenter.shinyapps.io/eDNAssay), a user-friendly tool for predicting qPCR cross-amplification.

Wolf Dispersal Patterns in the Italian Alps and Implications for Wildlife Diseases Spreading.

Wildlife dispersal directly influences population expansion patterns, and may have indirect effects on the spread of wildlife diseases. Despite its importance to conservation, little is known about dispersal for several species. Dispersal processes in expanding wolf (Canis lupus) populations in Europe is not well documented. Documenting the natural dispersal pattern of the expanding wolf population in the Alps might help understanding the overall population dynamics and identifying diseases that might be connected with the process. We documented 55 natural dispersal events of the expanding Italian wolf alpine population over a 20-year period through the use of non-invasive genetic sampling. We examined a 16-locus microsatellite DNA dataset of 2857 wolf samples mainly collected in the Western Alps. From this, we identified 915 individuals, recaptured 387 (42.3%) of individuals, documenting 55 dispersal events. On average, the minimum straight dispersal distance was 65.8 km (±67.7 km), from 7.7 km to 517.2 km. We discussed the potential implications for maintaining genetic diversity of the population and for wildlife diseases spreading.

Do metapopulations and management matter for relict headwater bull trout populations in a warming climate?

Mountain headwater streams have emerged as important climate refuges for native cold-water species due to their slow climate velocities and extreme physical conditions that inhibit non-native invasions. Species persisting in refuges often do so as fragmented, relict populations from broader historical distributions that are subject to ongoing habitat reductions and increasing isolation as climate change progresses. Key for conservation planning is determining where remaining populations will persist and how habitat restoration strategies can improve biological resilience to enhance the long-term prospects for species of concern. Studying bull trout, a headwater species in the northwestern USA, we developed habitat occupancy models using a data set of population occurrence in 991 natal habitat patches with a suite of novel geospatial covariates derived from high-resolution hydroclimatic scenarios and other sources representing watershed and instream habitat conditions, patch geometry, disturbance, and biological interactions. The best model correctly predicted bull trout occupancy status in 82.6% of the patches and included effects for: patch size estimated as habitat volume, extent of within-patch reaches <9°C mean August temperature, distance to nearest occupied patch, road density, invasive brook trout prevalence, patch slope, and frequency of high winter flows. The model was used to assess 16 scenarios of bull trout occurrence within the study streams that represented a range of restoration strategies under three climatic conditions (baseline, moderate change, and extreme change). Results suggested that regional improvements in bull trout status were difficult to achieve in realistic restoration strategies due to the pervasive nature of climate change and the limited extent of restoration actions given their high costs. However, occurrence probabilities in a subset of patches were highly responsive to restoration actions, suggesting that targeted investments to improve the resilience of some populations may be contextually beneficial. A possible strategy, therefore, is focusing effort on responsive populations near more robust population strongholds, thereby contributing to local enclaves where dispersal among populations further enhances resilience. Equally important, strongholds constituted a small numerical percentage of patches (5%-21%), yet encompassed the large majority of occupied habitat by volume (72%-89%) and their protection could have significant conservation benefits for bull trout.

Molecular species delimitation refines the taxonomy of native and nonnative physinine snails in North America.

Being able to associate an organism with a scientific name is fundamental to our understanding of its conservation status, ecology, and evolutionary history. Gastropods in the subfamily Physinae have been especially troublesome to identify because morphological variation can be unrelated to interspecific differences and there have been widespread introductions of an unknown number of species, which has led to a speculative taxonomy. To resolve uncertainty about species diversity in North America, we targeted an array of single-locus species delimitation methods at publically available specimens and new specimens collected from the Snake River basin, USA to generate species hypotheses, corroborated using nuclear analyses of the newly collected specimens. A total-evidence approach delineated 18 candidate species, revealing cryptic diversity within recognized taxa and a lack of support for other named taxa. Hypotheses regarding certain local endemics were confirmed, as were widespread introductions, including of an undescribed taxon likely belonging to a separate genus in southeastern Idaho for which the closest relatives are in southeast Asia. Overall, single-locus species delimitation was an effective first step toward understanding the diversity and distribution of species in Physinae and to guiding future investigation sampling and analyses of species hypotheses.

Occupancy modeling and resampling overcomes low test sensitivity to produce accurate SARS-CoV-2 prevalence estimates.

BACKGROUND: We evaluated whether occupancy modeling, an approach developed for detecting rare wildlife species, could overcome inherent accuracy limitations associated with rapid disease tests to generate fast, accurate, and affordable SARS-CoV-2 prevalence estimates. Occupancy modeling uses repeated sampling to estimate probability of false negative results, like those linked to rapid tests, for generating unbiased prevalence estimates. METHODS: We developed a simulation study to estimate SARS-CoV-2 prevalence using rapid, low-sensitivity, low-cost tests and slower, high-sensitivity, higher cost tests across a range of disease prevalence and sampling strategies. RESULTS: Occupancy modeling overcame the low sensitivity of rapid tests to generate prevalence estimates comparable to more accurate, slower tests. Moreover, minimal repeated sampling was required to offset low test sensitivity at low disease prevalence (0.1%), when rapid testing is most critical for informing disease management. CONCLUSIONS: Occupancy modeling enables the use of rapid tests to provide accurate, affordable, real-time estimates of the prevalence of emerging infectious diseases like SARS-CoV-2.

Identifying Candidate Genetic Markers of CDV Cross-Species Pathogenicity in African Lions.

Canine distemper virus (CDV) is a multi-host pathogen with variable clinical outcomes of infection across and within species. We used whole-genome sequencing (WGS) to search for viral markers correlated with clinical distemper in African lions. To identify candidate markers, we first documented single-nucleotide polymorphisms (SNPs) differentiating CDV strains associated with different clinical outcomes in lions in East Africa. We then conducted evolutionary analyses on WGS from all global CDV lineages to identify loci subject to selection. SNPs that both differentiated East African strains and were under selection were mapped to a phylogenetic tree representing global CDV diversity to assess if candidate markers correlated with documented outbreaks of clinical distemper in lions (n = 3). Of 54 SNPs differentiating East African strains, ten were under positive or episodic diversifying selection and 20 occurred in the clinical strain despite strong purifying selection at those loci. Candidate markers were in functional domains of the RNP complex (n = 19), the matrix protein (n = 4), on CDV glycoproteins (n = 5), and on the V protein (n = 1). We found mutations at two loci in common between sequences from three CDV outbreaks of clinical distemper in African lions; one in the signaling lymphocytic activation molecule receptor (SLAM)-binding region of the hemagglutinin protein and another in the catalytic center of phosphodiester bond formation on the large polymerase protein. These results suggest convergent evolution at these sites may have a functional role in clinical distemper outbreaks in African lions and uncover potential novel barriers to pathogenicity in this species.

Landscape genetics of wolverines (Gulo gulo): scale-dependent effects of bioclimatic, topographic, and anthropogenic variables.

Climate change can have particularly severe consequences for high-elevation species that are well-adapted to long-lasting snow conditions within their habitats. One such species is the wolverine, Gulo gulo, with several studies showing a strong, year-round association of the species with the area defined by persistent spring snow cover. This bioclimatic niche also predicts successful dispersal paths for wolverines in the contiguous United States, where the species shows low levels of genetic exchange and low effective population size. Here, we assess the influence of additional climatic, vegetative, topographic, and anthropogenic, variables on wolverine genetic structure in this region using a multivariate, multiscale, landscape genetic approach. This approach allows us to detect landscape-genetic relationships both due to typical, small-scale genetic exchange within habitat, as well as exceptional, long-distance dispersal among habitats. Results suggest that a combination of snow depth, terrain ruggedness, and housing density, best predict gene flow in wolverines, and that the relative importance of variables is scale-dependent. Environmental variables (i.e., isolation-by-resistance, IBR) were responsible for 79% of the explained variation at small scales (i.e., up to ~230 km), and 65% at broad scales (i.e., beyond ~420 km). In contrast, a null model based on only space (i.e., isolation-by-distance, IBD) accounted only for 17% and 11% of the variation at small and broad scales, respectively. Snow depth was the most important variable for predicting genetic structures overall, and at small scales, where it contributed 43% to the variance explained. At broad spatial scales, housing density and terrain ruggedness were most important with contributions to explained variation of 55% and 25%, respectively. While the small-scale analysis most likely captures gene flow within typical wolverine habitat complexes, the broad-scale analysis reflects long-distance dispersal across areas not typically inhabited by wolverines. These findings help to refine our understanding of the processes shaping wolverine genetic structure, which is important for maintaining and improving functional connectivity among remaining wolverine populations.

Pliocene-Early Pleistocene Geological Events Structure Pacific Martens (Martes caurina).

The complex topography, climate, and geological history of Western North America have shaped contemporary patterns of biodiversity and species distributions in the region. Pacific martens (Martes caurina) are distributed along the northern Pacific Coast of North America with disjunct populations found throughout the Northwestern Forested Mountains and Marine West Coast Forest ecoregions of the West Coast. Martes in this region have been classified into subspecies; however, the subspecific designation has been extensively debated. In this study, we use genomic data to delineate conservation units of Pacific marten in the Sierra-Cascade-Coastal montane belt in the western United States. We analyzed the mitochondrial genome for 94 individuals to evaluate the spatial distribution and divergence times of major lineages. We further genotyped 401 individuals at 13 microsatellite loci to investigate major patterns of population structure. Both nuclear and mitochondrial DNA suggest substantial genetic substructure concordant with historical subspecies designations. Our results revealed that the region contains 2 distinct mitochondrial lineages: a Cascades/Sierra lineage that diverged from the Cascades/coastal lineage 2.23 (1.48-3.14 mya), consistent with orogeny of the Cascade Mountain chain. Interestingly, Pacific Martes share phylogeographic patterns similar with other sympatric taxa, suggesting that the complex geological history has shaped the biota of this region. The information is critical for conservation and management efforts, and further investigation of adaptive diversity is warranted following appropriate revision of conservation management designations.

Integrative taxonomy refutes a species hypothesis: The asymmetric hybrid origin of Arsapnia arapahoe (Plecoptera, Capniidae).

Molecular tools are commonly directed at refining taxonomies and the species that constitute their fundamental units. This has been especially insightful for groups for which species hypotheses are ambiguous and have largely been based on morphological differences between certain life stages or sexes, and has added importance when taxa are a focus of conservation efforts. Here, we examine the taxonomic status of Arsapnia arapahoe, a winter stonefly in the family Capniidae that is a species of conservation concern because of its limited abundance and restricted range in northern Colorado, USA. Phylogenetic analyses of sequences of mitochondrial and nuclear genes of this and other capniid stoneflies from this region and elsewhere in western North America indicated extensive haplotype sharing, limited genetic differences, and a lack of reciprocal monophyly between A. arapahoe and the sympatric A. decepta, despite distinctive and consistent morphological differences in the sexual apparatus of males of both species. Analyses of autosomal and sex-linked single nucleotide polymorphisms detected using genotyping by sequencing indicated that all individuals of A. arapahoe consisted of F1 hybrids between female A. decepta and males of another sympatric stonefly, Capnia gracilaria. Rather than constitute a self-sustaining evolutionary lineage, A. arapahoe appears to represent the product of nonintrogressive hybridization in the limited area of syntopy between two widely distributed taxa. This offers a cautionary tale for taxonomists and conservation biologists working on the less-studied components of the global fauna.

An improved environmental DNA assay for bull trout (Salvelinus confluentus) based on the ribosomal internal transcribed spacer I.

The majority of environmental DNA (eDNA) assays for vertebrate species are based on commonly analyzed regions of the mitochondrial genome. However, the high degree of mitochondrial similarity between two species of charr (Salvelinus spp.), southern Dolly Varden (S. malma lordii) and bull trout (Salvelinus confluentus), precludes the development of a mitochondrial eDNA assay to distinguish them. Presented here is an eDNA assay to detect bull trout based on the first ribosomal internal transcribed spacer (ITSI), a nuclear marker. This assay successfully detects bull trout and avoids detecting Dolly Varden as well as brook trout (S. fontinalis), Arctic char (S. alpinus), and lake trout (S. namaycush). In addition, this assay was compared with an extensively used mitochondrial bull trout assay and it was found that the ITSI-based assay produced higher detectability. Our results suggest this assay should out-perform the published mtDNA assay across the range of bull trout, while the added specificity allows reliable bull trout detection in areas where bull trout co-occur with other charr such as Dolly Varden. While clearly a superior assay in this instance, basing assays on ITSI is not without problems. For vertebrates, there are fewer ITSI sequences available than commonly sequenced regions of the mitochondrial genome. Thus, the initial in silico screening of candidate assays must be preceded by much more extensive sampling and sequencing of sympatric or closely related taxa. Further, all copies of the internal transcribed spacers within an individual may not be identical, which can lead to complications. Lastly, the copy number for ITSI varies widely across taxa; the greater detectability associated with this assay cannot be assumed for other species.